Tuning the electron-phonon coupling in multilayer graphene with magnetic fields 
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Magneto Raman scattering study of the F<2 g optical phonons in multi-layer epitaxial graphene 
grown on a carbon face of SiC are presented. At 4.2K in magnetic field up to 33 T, we observe 
a series of well pronounced avoided crossings each time the optically active inter Landau level 
transition is tuned in resonance with the E2 9 phonon excitation (at 196 meV). The width of the 
phonon Raman scattering response also shows pronounced variations and is enhanced in conditions 
of resonance. The experimental results are well reproduced by a model that gives directly the 
strength of the electron-phonon interaction. 

PACS numbers: 73.22.Lp, 63.20.Kr, 78.30.Na, 78.67.-n 
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Lattice vibrations in solids can be effectively modified 
via their coupling to electronic excitations, as, for ex- 
ample, evidenced by observations of Kohn anomalies in 
metals 0, Hj] , of coupled phonon-plasmon modes in polar 
semiconductors 0, |j] or of different phonon spectra in 
metallic and semiconducting carbon-nanotubes The 
electron-phonon interaction is currently intensively stud- 
ied in graphene @, 0, B @| which is a two-dimensional 
crystal of carbon atoms arranged in a honeycomb lattice 
and a semimetal with characteristic dispersions of elec- 
tronic states displaying Dirac cones near the Fermi en- 
ergy 

num. The 

case of the long wavelength optical E 2g - 
phonons at the T-point of the Brillouin zone, which corre- 
sponds to the relative displacement of two nonequivalent 
carbon atoms in the unit cell of graphene @, 0] , is of par- 
ticular interest. The perturbation due to this displace- 
ment is effective in inducing the direct (Ak = 0) elec- 
tronic transitions across the Dirac point: E2 S -phonons 
efficiently couple to those low energy interband excita- 
tions (US] that are unique to graphene. The spectrum 
of these excitations (its 2Ep low energy onset) can be 
modified by tuning the Fermi energy Ep. This was 
achieved in gated graphene flakes on Si/SiC>2 substrates 
where the electrically modified E2 9 -phonon spectrum was 
traced with Raman scattering methods [1, 0| . 

The spectrum of the graphene E2 9 -phonon is expected 
to be even more severely modified by applying a magnetic 



field perpendicular to the 2D plane [12lll3l|. i.e., when a 
continuous spectrum of electronic excitations is trans- 
formed into a series of quasi-discrete inter Landau level 
excitations characteristic of a 2D system. In conditions of 
Landau quantization, the electron phonon coupling has 
a resonant character which is expected between the E2 ff - 
phonon and properly selected inter-Landau level excita- 
tions. The observation of the effects of magneto-resonant 
electron-phonon coupling in graphene structures, which 



we present in this letter, has been an experimental chal- 
lenge 12711 aiming at verification of theoretical predic- 
tions [l2l . [l3| and eventual confirmation of the conclu- 
sions drawn from zero- field measurements [H, Q . 

We report here on magneto-Raman scattering studies 
of the E2 g -phonon band of multilayer epitaxial graphene 
on the carbon face of a SiC substrate (MEG), in fields 
up to 33 T, and low, liquid helium temperatures. The 
measured E2 ff -band is observed to be composed of two 
components: one fixed and another oscillating with the 
applied field. The extracted oscillatory component of the 
E2 9 -band shows the theoretically anticipated magneto- 
phonon multi resonance [l2l [l3[ . Pronounced variations 
of the energy position and width of this line can be un- 
derstood in terms of a series of avoided crossings each 
time the E2 9 -phonon energy tends to match the energy 
of the specific (defined by appropriate selection rules) in- 
ter Landau level excitations. Comparison of the experi- 
mental results to existing models allows us to analyze the 
characteristic strength of the electron-phonon interaction 
in the investigated graphene system. Possible origin of 
the field-independent component of the E2 9 -band is dis- 
cussed. 

Raman scattering spectra were measured using the 
Ti:Saphire laser setup, tuned at accurately controlled 
wavelength in the range ~ A=720 nm, in order to mini- 
mize the superfluous Raman signal of optical fibers which 
were used to transmit the excitation light to, and to col- 
lect the Raman photons from the sample. The sample 
was immersed in a helium gas kept at T= 4.2 K and 
placed in a resistive magnet delivering fields up to 33 T. 
The non-polarized Raman scattering spectra were mea- 
sured in nearly back-scattering Faraday geometry: the 
arrangement of the excitation and collection fibers (both 
with core diameters of 200 /im) with respect to sample 
surface is sketched in the inset of Fig.QJt. Typical excita- 
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tion power was ^100 mW, spread over 600 /im -diameter 
laser spot on the sample. The collected light was dis- 
persed with a single grating spectrometer (spectral reso- 
lution AA=0.3nm) equipped with nitrogen cooled CCD 
detector and band pass filters were used to reject the 
stray light. 

The MEG structures which we have used for experi- 
ments were fabricated by thermal decomposition of the 
carbon face of a 4H-SiC substrate Two investi- 

gated structures, both with a relatively large number of 
graphene layers, show similar data and we present here 
the results obtained for the sample which contains ~ 70 
layers. We used strongly graphitized samples to enhance 
the Raman scattering signal, particularly because the 
E2 9 -phonon band in MEG samples appears on the back- 
ground of the two-phonon Raman scattering processes in 
the SiC substrate. A number of exp erimental studies, 
including magneto-absorption ID, 16| and micro-Raman 
scattering measurements 17j, show that Dirac-like elec- 
tronic bands, which are genuine of a single graphene 
sheet, persist in MEG structures and that a large part of 
these layers is practically neutral, with carrier densities as 
low as 5 x 10 9 cm -2 and fairly high mobilites up to 250000 
cm 2 /(V.s) [la]. The electronic bands in MEG grown on 
the carbon face of SiC are alike those of a single layer 
because of preferentially rotational and not Bernal-type 
layer stacking in this material HI, ljj, 2(], 21 1. Although 
we believe that the majority of the structure probed in 
our experiments displays neutral graphene-like electronic 
system, the MEG samples and in particular the strongly 
graphitized specimens are not expected to be uniform on 
the scale of 600 /im defined by the actual size of the laser 
spot on the sample. For example, the presence of Bernal 
stacked minority inclusions in these structures has been 
seen with micro- Raman scattering experiments Uj . 

The representative Raman scattering bands observed 
for our ~ 70 layers MEG sample in the absence of mag- 
netic field are shown in Fig. QJi. The so-called "2D" 
band, which is frequently analyzed to characterize differ- 
ent graphene structures [22|, H^], appears at 2646 cm -1 
as a 35 cm _1 -wide, slightly asymmetric line. A pure 
single Lorentzian shape of the "2D" line is a signature 
of the electronic system with simple Dirac-like electronic 
states, such as those characteristic of graphene or of the 
majority part of the MEG structures. We understand 
that the observation of the high energy asymmetry in the 
" 2D" line under the present experimental conditions is a 
consequence of signal averaging from a relatively large 
portion of the sample, which includes other than simple 
graphene-like structures (i.e., Bernal stacked residues). 

The focus of this work is on the E2 S -phonon line which 
clearly dominates over the background of two-phonon sig- 
nal from SiC substrate and peaks at 1586.5 cm -1 in the 
Raman scattering spectrum measured at zero magnetic 
field (Fig. QJ,). As shown in Fig. QJ> with background 
corrected spectra, the E2 9 -phonon line is visibly affected 
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FIG. 1: (a) The E2 9 and 2D Raman scattering bands mea- 
sured at T= 4.2 K and A = 720.7 nm laser excitation at zero 
magnetic field. The scheme of the experimental configura- 
tion is shown in the inset, (b) Background corrected spectra 
of the E2 9 band, measured at different magnetic fields. The 
solid red lines represent the fitted two-Lorentzian component 
shape. One of this component is field independent and fixed 
to solely reproduce the spectrum at B = 26 T. The oscillatory 
with field component is shown with dotted black lines. 



by the application of the magnetic field. A more de- 
tailed data inspection shows however that the E2 9 line 
is practically field independent at high field region above 
B> 24 T. Hence, one may anticipate that the observed 
magneto-Raman spectrum of the E2 ff -phonon is com- 
posed of two components: one oscillating with field and 
the second one which is field independent and dominant 
at B > 24 T. Indeed, as illustrated in Fig. [lb, the E 2ff 
spectrum measured at any magnetic field can be fairly 
well fit with two-Lorentzian functions (solid red lines). 
The field independent component I26 (hu) has been set to 
reproduce the spectrum at B = 26 T. The resulting oscil- 
latory component is shown with dashed lines in Fig. [TJd. 
To present the full data collection, the I26 (hi/) Lorentzian 
function has been extracted from each measured spec- 
trum and the resulting differential spectra (oscillatory 
part) are shown in Fig. [2l in form of the intensity false 
color map. Full points in this figure denote the peak po- 
sition (center) of the Lorentzian functions which fit the 
oscillatory component. 

The magnetic field evolution of the spectra shown in 
Fig. [5] can be clearly seen as resulting from a series of 
avoided crossings between the E2 9 phonon excitation and 
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FIG. 2: a) Color map of the magneto-oscillatory component 
of Raman scattering spectra of E2 9 band phonons as a func- 
tion of the magnetic field measured at T = 4.2 K under 
A=720.7 nm laser excitation. The extracted peak position of 
this line is shown with full dots. Their size is proportional to 
the line amplitude. Solid lines Tfc represent the energies of the 
series of inter Landau level transitions:!/. 



fc.(-fc-i) 
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which couple to E2 S phonon (e = 1.02 x 10 6 m.s 1 is assumed), 
b) Zoom on the to 10 T range of magnetic fields. 



the specific inter Landau level electronic transitions. The 
energy ladder of Landau levels (L n ; n = 0, ±1, ±2, ...) of 
the two-dimensional Dirac-like electronic system is given 
by E„ = sgn(n)cy/2ehB\n\, where the measured slope of 
the Dirac cone (Fermi velocity) c = 1.02 • 10 6 m/s in the 



investigated structures 



velocity) 

HEfl. As illustrated in Fig. El 
the E 2g -phonon couples to a series of L n — > L rn tran- 
sitions which fulfill the \n\ — \m\ = ±1 selection rules 
and therefore appear at energies = («Jk + \Jk + X)E\ 
(k = 0,1,2,..). Notable, the same selection rules apply 
to optically active inter Landau level transitions which 
are observed in far-infrared magneto-transmission spec- 
tra [H, H 0, El- The amplitude of the observed "an- 



ticrossings" grows with the magnetic field because, as a 
rule, the splitting (A) between coupled modes is propor- 
tional to the square root of the oscillator strength of the 
uncoupled excitations and to the square root of interac- 
tion parameter A. In our case, these are electronic ex- 



citations which become enhanced by the magnetic field. 
This is primarly due to the increasing with B-field Lan- 
dau level occupation [eB/h degeneracy). If broadening is 
neglected we roughly expect A ~ \J \B and derive more 
rigorously A = \f2\E\, for our neutral Dirac-like 
system. 

For the sake of more quantitative data analysis we con- 
sider both the characteristic phonon energy as well as its 
broadening which is another inherent element of mode 
coupling effects. In Fig. Owe re-plot the center and add 
the extracted linewidth (HWHM) of the Lorentzian peak 
which fits the magneto-oscillatory component of the mea- 
sured E 2g -line. The magneto-oscillations in the peak po- 
sition are accompanied by oscillations in the linewidth 
which, as expected, is enhanced under resonant condi- 
tions. To reproduce the data of Fig. [3l we adopt the 
approach of Ando [lj| to our neutral Dirac-like system 
and derive the phonon energy e and broadening param- 
eter F by extracting e = e — iT from the equation which 
defines the poles of the phonon Green's function: 
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where eo stands for the phonon energy of the neu- 
tral system at B=0 T and 8 accounts for the broaden- 
ing characteristic for electronic excitations. The mea- 
sured linewidth has been assumed as a convolution sum 
\/5q + T 2 , where 5q accounts for other, than electron- 
phonon coupling, broadening mechanisms. 

To calculate the energy and broadening of the E 2g — 
phonon spectrum we fix the value of c = 1.02 x 10 6 m/s 
as derived from far-infrared magneto-absorption experi- 
ments on MEG samples 15, III ■ Working effectively with 
only two adjustable parameters S and A (Sq = 3.7 cm -1 is 
the background correction to the linewidth oscillations) 
the experimental data can be fairly well reproduced. Op- 
timal choice of parameters is: <5 = 90 cm^ 1 and A = 
4.5 x 10~ 3 . The extracted strength of electron-phonon 
coupling is in excellent agreement with estimations of A in 
the range of 4—5.5 x 10 -3 from the detailed analysis of the 
dispersions of the E 2s in graphene (25| and/or measure- 
ments of electrically tuned E 2g -phonon in graphene 0, 0] . 
We also note that the width S which is measured in be- 
tween well separated resonances, i.e., when the effects of 
electron-phonon interaction are switched off, appears to 
be practically identical to its value measured under cor- 
responding conditions (2Ef > eo) on graphene flakes on 
Si/Si0 2 [8]. The parameter S accounts for broadening of 
electronic states and its actual value reasonable matches 
the width of magneto-absorption transitions expected in 
our sample at high magnetic fields [26j]. For simplicity 
6 is assumed constant, thought its increase with field 
(~ V~B) which is expected from magneto-absorption ex- 
periments [l6| could improve the data simulation (larger 
broadening of the high field resonance). 
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FIG. 3: Magnetic field evolution of the peak position (upper 
part) and the linewidth (HWHM) of the magneto-oscillatory 
component of the measured E2 9 Raman band. Solid lines rep- 
resent the result of modelling of the data (with the parameters 
specified in the figure) along the procedure described in the 
text. 



The puzzle of the presented results is the origin of the 
field independent component of the E2 9 spectrum which, 
we note, represents the significant portion of the total 
spectral weight. The Bernal-stacked inclusions are first 
possible candidates to account for this contribution (we 
have checked that in fields up to 14T the magneto-Raman 
spectrum of the F,2 g band of natural graphite sample is 
practically field independent). However, inspecting the 
shape of the measured 2D band, we can rule out that the 
contribution from these inclusions is as large as 65% (ex- 
tracted for the B = 0T spectrum) . Another candidate is 
a signal from highly doped layers close to the interface. 
This is also scarcely probable since such signal should be 
significantly shifted in frequency from that of the major- 
ity neutral graphene layers, which is not seen. We specu- 
late the field independent component of the E2 ff line is an 
inherent property of the graphene-like systems [27] what 
remains to be understood. 

In conclusion, we have investigated the magnetic field 
evolution of the Raman scattering response of the F,2 g 
phonons in a graphene-like system in fields up to 33 T. A 
pronounced avoided crossing behavior of the phonon en- 
ergy and oscillations of the line width are observed each 
time the optically active inter Landau level excitation is 
tuned in resonance with the phonon energy. The exper- 
imental results are well accounted for by existing the- 



oretical models 12j, [13J of electron-phonon coupling in 
graphene, with the estimated electron-phonon coupling 
constant A = 4.5 x 10 -3 . 

We would like to acknowledge fruitful discussions with 
V. Fal'ko and M.O. Goerbig. We acknowledge T. Kaz- 
imierczuk for help with the data analysis and A. Wys- 
molek for his help in the initial stage of this experiment. 
Part of this work has been supported by ANR projects 
ANR-08-JCJC-0034-01, ANR-06-NANO-019 and by the 
French Embassy (Partner University Fund). One of us 
(P.K.) is financially supported by the EU under FP7, con- 
tract no. 221515 'MOCNA'. M. A. thanks the Cariplo 
Foundation (project QUANTDEV), MICINN (Project 
MOSAICO) and JCYL SA052A07 for support. 



[1 

[2. 
[3 

[4. 
[5 
[6] 
[7] 

[8 

[9. 
[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 
[18 

[19 

[20 

[21 
[22 
[23 
[24 

[25 

[26 
[27 



197004 (2004). 
Rev. 16, 999 



W. Kohn, Phys. Rev. 2, 393 (1959). 
A. Q. R. Baron et al, Phys. Rev. Lett. 92, 
A. Mooradian and G. B. Wright, Phys. 
(1966). 

A. Wysmolek et al, Phys. Rev. B 74, 165206 (2006). 

S. Reich et al, Carbon Nanotubes (2004). 

T. Ando, J. Phys. Soc.Jpn. 75, 124701 (2006). 

A. H. C. Neto and F. Guinea, Phys. Rev. B 75, 045404 

(2007). 

J. Yan et al, Phys. Rev. Lett. 98, 166802 (2007). 

S. Pisana et al., Nat. Mat. 6, 198 (2007). 

A. K. Geim and K. S. Novoselov, Nature Materials 

(2007). 

A. H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 
Ando, J. Phys. Soc.Jpn. 76, 024712 (2007). 
O. Goerbig et al, Phys. Rev. Lett. 99, 087402 (2007). 
Berger et al., J. Phys. Chem. 108, 19912 (2004). 
L. Sadowski et al, Phys. Rev. Lett. 97, 266405 (2006). 
Orlita et al., Phys. Rev. Lett. 101, 267601 (2008). 
Faugeras et al., Appl. Phys. Lett. 92, 011914 (2008). 
J. Hass et al, Phys. Rev. Lett. 100, 125504 (2008). 
J. L. dos Santos, N. M. R. Peres, and A. H. C. Neto, 
Phys. Rev. Lett. 99, 256802 (2007). 

S. Latil, V. Meunier, and L. Henrard, Phys. Rev. B 76, 
201402(R) (2007). 

L. Magaud et al, Phys. Rev. B 79, 161405(R) (2009). 

A. C. Ferrari et al, Phys. Rev. Lett. 97, 187401 (2006). 

D. Graf et al, Nano Lett. 7, 238 (2007). 

M. L. Sadowski et al., Int. J. of Mod. Phys. B 21, 1145 

(2007). 

M. Lazzeri and F. Mauri, Phys. Rev. Lett. 97, 266407 
(2006). 

P. Plochocka et al, Phys. Rev. Lett. 100, 087401 (2008). 
The experimental approaches to trace the magneto- 
phonon resonance effects in graphene flakes on Si/Si02 
substrates have been presented by J. Yan et al. at the 
APS 2009 March meeting in Pittsburg and by S. Remi et 
al. at the KIAS Graphene meeting in Seoul. Notably, S. 
Remi et al. has reported the splitting and not a simple 
shift of the E2 9 phonon in the Raman spectrum at B = 
11 T of likely quasi-neutral graphene. 



